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Abstract 

Using the Bohzmann equation with a Langevin-hke term describing the 
stochastic force in a baryon-photon plasma, we investigate the influence of 
the incoherent electron-photon scattering on the subhorizon evolution of the 
cosmic microwave radiation. The stochastic fluctuation caused by each col- 
lision on average is found to be small. Nevertheless, it leads to a significant 
Brownian drifting of the phase in the acoustic oscillation, and the coherent 
oscillations cannot be maintained during their dynamical evolution. As a 

consequence, the proposed Doppler peaks probably do not exist. 
PACS numbers: 98.70.Vc, 98.80.Es, 05.40.+j 
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The anisotropies of the cosmic microwave background (CMB) provide a key to the un- 
derstanding of the origin of primordial fluctuations and the thermal history of the early 
universe. The spectrum of the CMB anisotropy on large angular scales has been found to 
be consistent with the inflationary scenario of the early universe [|I|. It is generally believed 
that many cosmological parameters can be determined from the fine structures in the power 
spectrum of CMB anisotropy. Among them, the possible existence of Doppler peaks - the 
peaks in the CMB anisotropy spectrum on an angular scale of about one degree or less has 
attracted much attention 0. 

It has been shown theoretically that the amplitude and the position of these Doppler 
peaks are functions of the spatial curvatures, mass density, reionization time, cosmological 
constant etc. Thus, a precision measurement of the Doppler peaks may provide us an 
effective tool to determine various cosmic parameters. However, the observed amplitudes of 
the CMB anisotropy on degree-scale have not so far been very conclusive in determining the 
existence of the Doppler peaks. Indeed, some observations seems to exhibit high amplitudes 
as expected in a Doppler-peak scenario, while others show no peak amplitudes One may 
assert at this stage that the expected peaks have not yet been clearly identified in current 
data, but would be determined by a new generations of the CMB anisotropy observations. 

In this letter, we shall take a different approach to reexamine the theoretical foundation 
of the prediction of the Doppler peaks. In particular, we shall argue based on a stochastic 
Boltzmann equation that the coherence of acoustic oscillations in the baryon-photon plasma 
will be largely disturbed, and may even be totally erased if the stochastic force of the 
incoherent scattering is included. In the standard theory of the CMB evolution 0, the 
acoustic oscillations in the baryon-photon plasma on subhorizon scales are treated coherently, 
i.e. different modes are frozen at different phases of their oscillation. The position of 
the Doppler peaks is then calculated by the phase at the recombination. However, the 
inclusion of correlation due to the stochastic term in the baryon-photon plasma will lead to 
a phase randomization. The coherence of the oscillations could be maintained if there is a 
mechanism for providing a negative entropy current to prevent the decoherence due to the 
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phase randomization (like a laser). Unfortunately, no such mechanism exists in the epoch 
of recombination, while the existence of a stochastic force term in the kinetic equation 
is inevitable in a system with dissipation. We shall show that the phase randomizations 
introduced by the incoherent electron-photon scattering are, indeed, substantial, and as a 
result, the Doppler peaks will be erased. Instead, one expects a large dispersion of the 
amplitudes in the CMB power spectrum due to different realizations on subhorizon scales, 
reflecting the stochastic nature of the phases. We suggest that the dispersed observations of 
the CMB angular power spectrum on degree-scale, though very coarse at this stage, are in 
support of our view. 

Let us use conventional notations in the theory of the CMB anisotropy [Q. Since the 
coherent oscillations of the Doppler peaks act only on subhorizon scales, the choice of gauge 
is irrelevant for our calculation. The evolution of the photon distribution function f{t, x, p) 
is calculated by the Boltzmann equation 

dt~ dt dt dx^ dt dpo dt d-f' ^ ' 

where 7^ is the direction cosines of pi with respect to the corresponding spatial coordinate. 
The left hand side of eq.(l) describes the free-streaming, and the right hand side is the 
collision term given by 

C[f]= y'rfqf/q'rfp'H^(p,q,p',q') 

x{(7(t,x,q')/(t,x,p')[l + /(t,x,p)] 

-^(t,x,q)/(t,x,p)[l + /(t,x,p')]} (2) 

where g is the electron distribution function, and the collision rate W is determined by the 
Compton scattering between electron and photon from state (q', p') to (q, p). 

It is well known that the Boltzmann equation in (1) is derived under the assumption of 
a molecular chaos and is applicable if the fluctuations caused by the incoherent collisions 
are negligible. These fluctuations give rise to a stochastic force term in hydrodynamics [Q, 
governed by the fluctuation-dissipation theorem. Similarly, these fluctuations can be taken 
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into account by an additional Langevin-like force term in the Boltzmann equation [Q. Eq.(l) 
should then be replaced by 

| = q/]+r(t,x,p) (3) 
where the stochastic force, r, is characterized by its correlations: 

(r(t,x,p)) = 0, (4) 

(r(ti,Xi,pi)r(t2,X2,P2)) = ^N6{ti - ta) 
x5(xi-X2) J dqdq^dpdp'W{q,p,c^, p')A(pi) 

X A(p2)(7(ti, xi, q)/(ti, xi, p) [1 + /(ti, xi, p')] (5) 

where (...) is an average over the stochastic effect (i.e. over different realizations), is the 
total number of photons, and the function A(P) is defined by 

A(P) = 5(P - q) + (5(P - p) 

-^(P-q')-^(P-p')- (6) 

It is important to note that eqs.(3)-(6) are applicable not only in a linear approximation but 
also in the nonlinear region |^,^. 

In principle, to study the fluctuations in the baryon-photon plasma, we should also 
consider the stochastic terms in the equation of baryonic matter. However, since the linear 
fluctuations given by independent Gaussian random "forces" are additive, the stochastic 
terms in the baryonic equation will increase the effects of fluctuations considered in this 
paper. To illustrate the main effect of electron-photon scattering in the presence of the 
Doppler peaks, we shall only consider the contribution from stochastic term r in photon's 
equation. In this case, the electron distribution (^(t, x, q) can be treated as an external 
source. For a Thomson scattering, photons do not exchange energies with electrons. The 
q-distribution of electrons do not involved in the evolution, and we need only the number 
density distribution of electrons Ue = J (iq5'(t, x, q). In this case, the distribution of the 
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photon energy p is also unchanged, the perturbation of the CMB can be described by the 
anisotropy of the brightness temperature 

where p-y = (7r^/15)T^ is the mean energy density of photons. In terms of 6, the photon 
distribution function can be approximately expressed as 

/(i,x,p) = /T(p)[l + 4e(t,x,7)] (8) 

where fx = l/[exp(p/T) — 1]. 

From eqs.(l) and (8), one finds that © should satisfy 

|_(e + *) + y|-(e + *) + ,.|-e = 

T(e„-e+7ii.* + ^Qe) + B (9) 

where $ and are the Newtonian potential and the space curvature perturbation, respec- 
tively. Tj = /(I + z)dt is the conformal time, z the redshift, vi, the baryonic (fluid) velocity, 
and r = ngUT the optical depth, err being the Thomson cross section. Q is the projection 
operator of quadrupole, defined as Q(7,7') = Yl^=-2^2m{l){^l^'^) I ^^'^2m(7')- Obviously 
QQ=Q. The isotropic component 6o(r/,x) is given by ©o — O©, where O = (l/47r) / dQ is 
the monopole projection operator, and 00=0. The term ©o appearing on the right-hand 
side of eq.(9) indicates that without an external driving force (such as Vb) the isotropic state 
is the "equihbrium" state in the kinetic evolution of the Thomson scattering. 

The stochastic term in eq.(9) is i? = (l/47r^p^) / rp^dp, and its correlation function is 
given by 

(i?(r7i,xi,7i)i?(r/2,X2,72)) = 
3 N 



87r^ lOTT^P^ 

x[^47r(5(7i-72)-^(l + cos2/3)] 



X J p'^MfiVi,^!,?, 7i) 

+2/(771, Xi,p,7i)/(772,X2,P, 72) +/(?72,X2,P, 72)] (10) 



where (3 is the angle between 71 and 72. Because eq.(9) is hnear in ©, it reduces to eq.(l) by 
taking an average over the stochastic effect. Thus, the calculation based on eq.(l) is actually 
only for (0) but not for 0, i.e. the fluctuations due to the stochastic terms are entirely 
overlooked in eq.(l). In the incoherent processes, the linear fluctuations, 5Q — Q — (0), 
caused by stochastic force R may not in general play a very important role. However, these 
fluctuations lead to "forgetting history" . Thus, the existence and maintenance of coherence 
should be seriously reconsidered. 

Let us calculate the linear fluctuations, ^0(77, x, 7), at a given time 77 around a solution 
to eq.(l). It is governed by 

id , d ■. d 

{ h 7 r + 7* r 

[ dr) dx^ dY 

+T[l + ^Q]}(50-<^0o)=i?' (11) 

where R' — {1 — 0)R. The projection operator arises from the isotropic term 0o in eq.(9). 
The solution to eq.(12) can be generally expressed as 



50(77, X, 7) - SOoir), X, 7) = / dXe' 

Jo 

xR'(r) -X,x'- I dX'Y, y - r dAV) (12) 
Jo Jo 



-Ar[l+(1/10)Q] 

rX 



Obviously, both sides of eq.(13) become zero upon an integration (l/47r) / dfl because 0(1 — 
O) =0. 

The correlation functions of fluctuations, such as (50(771, xi, 71)50(772, X2, 72)), can be 
calculated from eqs.(lO) and (12). For a flat universe, the light path is a straight line in 
comoving coordinates x, and thus 7* = 0. The fluctuation for the mode (i,k) with / > is 
then 

(50,(77, k)50;(77,k)) = 
— J (ixi(ix2 exp{-ik ■ (xi - X2)} 

dQidQsPKko • li)Pii'ko ■ 72) 



(2/ + 1) 



647r 



J dX,dX2exp{-TX,[l + Qi] - tA2[1 + Q2]} 



{R'{r] - Ai, xi - Ai7i, -fi)R'{r] - A2, Xa - A272, 72)) 

where ko = ^ volume being considered. Since (r/k) » 1, the above equation 

can be simphfied by using eqs.(8) and (10) as 

0.531 X 675 
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x^^^^ / d^^dn^Piiko ■ 7i)^/(ko • 72) 
J TrfAe-'^^(l-0i)(l-02) 
x[l + (e^^/^° - l)Qi][l + (e^^/^° - 1)Q2] 
x[^47r5(7i-72)-^(l + cos2/3)] 
xi| cix(e(r7-A,x,7i)e*(77-A,x,72)) (13) 

where we used TV = {2C{3)/7r^)VT\ J p'^dpf^ = 2[C(2) - C(3)]r^ = 0.844r5, and e"^ - 
[1 + (e" — 1)Q]. The monopole projection factors (1 — Oi) and (1 — O2) can be replaced 
by a unity operator when I > 0. The two quadrupolc projection terms Qi and Q2 can be 
neglected because the term e~^^'^ picks up the main contribution in the integral from the 
region At < 1. Using the following mode decomposition 



V 
2^ 



i I dx(e(77,x,7i)e*(77,x,72)) = 

roo 1 

I k'dkY:^^{\e,{v,k)\')Pi'hi-^2) (14) 



we have finally 



(|5eKr/,k)H = 3.91x10-^ 
-| E [5.,o+|(/ + i^,0,2,0|/ + i/,0)|2] 



X 



where the overline above {\Qi{ri, k)\)^ denotes an average over the region from 77 — (1/r) 
to 77, and {I + v, 0, 2, 0\l + u, 0) are the Clebsch-Gordon coefficients. For large I, the second 



(negative) term in the bracket eq.(16) is completely negligible. The first term is independent 
of I and k. This is expected since the correlation function of the corresponding stochastic 
force is isotropic in 7-space and uniform in fc-space. 
Using the expression for a temperature perturbation 

'AT\2 V r 1 



one has 



1.30 X 10"^|AT/T|^ (17) 



21 + 1 

The I summation in eq.(16) should run up to Imax where the electron-photon collision is 
frequent enough, i.e. to the scale about the photon mean free path 1/r. Therefore, AT/T 
in eqs.(16) and (17) should not be confused with {AT/T)obs given by observations with 
low resolutions, i.e. their window functions are on scales much larger than Imax- For the 
observation of CMR anisotropy with a resolution /, the fluctuation is 
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\T') t5o2/' + i; 

AT 

~ 3.61 X IQ-^Vl—. (18) 

Therefore, the fluctuations caused by the stochastic force R are generally small, except for 
the case of a very high resolution observation. 

However, the average of the fluctuations in eqs.(15) and (18) is essentially over one 
collision time 1/r. In relation to coherent processes, we should study the cumulative effect 
of the fluctuations over the entire period during which the coherence is to be maintained [0 . 
For the Doppler peaks, the period is from the time when the perturbation enters horizon to 
the time of recombination, i.e. the duration of the subhorizon evolution before recombination. 
The cumulative effect can be easily described by the phase fluctuations of the oscillation. 
Use the expression Qi{ri,k) = \Qi\e^'^', where |9;(?7, fc)| and (piijj.k) are the amplitude and 
phase, respectively. The equations of 6/(?7, k) can then be derived from eq.(9) as 
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eo = -^ei-$ + i?o (19) 

01 = (^Go + ^ - ^62) - r(0i - H) + Ri (20) 

-re, + i?/ (/>2) (21) 



where 



J dQPi{ko ■ 7) J d^e-'^ ^'R (22) 



V Stt 

If the hierarchy of (21)-(23) is cut off at 2/-th order, they correspond to the equations for a 
system consisting of / coupled oscillators. Without the stochastic terms Ri, the oscillations 
are coherent and the phases of the oscillations are completely fixed by initial conditions. 
The terms Ri lead to a phase randomization. 

We shall first calculate the phase fluctuations raised by Rq. As the coupling between 
the photon and the baryon is tight and the peculiar gravitational potential is weakly time- 
dependent in matter-dominated regime 0, the evolution of 60 is dominated by the phase 
evolution (polrj), which is approximately described by a WKB-like equation as 

dr^ - ' 2(60) kc, dri ' ^ ^ 

where Cg denotes the sound speed of the baryon-photon plasma, which is ~ before 
recombination. The stochasticly averaged solution of eq.(23) is 0o(^) = kcgdrj' + (pQ^rjen) , 
where rjen denotes the time when the mode k enters the horizon, and (poiVen) is the initial 
phase. The position of the Doppler peaks for adiabatic perturbations is approximately 
determined by a phase relation 0o('7re) = nir, where n is integer. Eq.(23) shows also that 
the phase fluctuation due to Rq in the period from rjen to rjre is 



r,.„ '2{eo)kcs dri' ' ^ ^ 



Generally, the term Ri leads to a phase fluctuation at least of 



t)i ~ 



r,e„ '2{Qi)kcs dr]' ' ^ ^ 
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Considering that the derivative {l/kcs)d/dr] in eq.(25) contributes a factor of order 1, and 
using eqs.(lO), (15) and (17), one obtains the mean phase fluctuation as 



The factor {AT/Tf/{\5ei{r],k)\'^) is no less than 1, and therefore, the RMS of the phase 
fluctuation can be estimated as 



where N^. = T{rjre — Ven) is the mean number of collisions in the entire period of subhorizon 
evolutions of mode k before the recombination. Therefore, the behavior of the stochastic 
fluctuation of the phase for the Z-oscillation is just like a Brownian drifting: the number of 
collisions corresponds to the number of steps in the random walk, and the mean shift per step 
is about 1.8x 10"^"\/2/ + 1. For Doppler peaks, we have / > 100, and N > rlrjre-rjen) > 10^. 
Therefore, we conclude that the RMS of the Brownian phase drifting due to Thomson 
scattering is order 1, significant enough to disturb the coherence. 

In addition to the Thomson scattering, there are other stochastic forces in the baryon- 
photon plasma, such as non- Thomson terms of the Compton scattering, the stochastic terms 
in the hydrodynamical equation of baryons. Fluctuations from different sources mostly are 
additive. Our estimation on the phases drifting is thus a conservative one. The proposed 
Doppler peaks in the CMB angular power spectrum probably do not exist. Different causal 
areas at the recombination can be considered as independent realizations of the stochastic 
force. The dispersion of the currently observed CMB anisotropics on angular scale of one 
degree is consistent with the scenario of the Brownian drift of the phases. 
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